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The gauge-invariant, nonlocal, dynamical quark model is proved to generate the minimal vertices 
which satisiy the Ward-Takahashi identities. In the chiral limit, the momentum-dependent quark 
self-energy results in a flat-like form with some end point 5— funtions for the light-cone pion dis- 
tribution amplitudes, similarly found in the Nambu Jona-Lasino model with constant constituent 
mass. The leading order nonzero pion and current quark masses corrections lead concave type 
asymptotic-like form modifications to twist-2 pion distribution amplitude with end point pillars and 
twist-3 tensor pion distribution amplitude above the fiat-like form backgrounds. A by-product of 
our investigation shows that the variable u appearing in pion light-cone distribution amplitudes is 
just the standard Feynman parameter in the Feynman parameter integrals; also chiral perturbation 
works well for these amplitudes. 

PACS numbers: ll.lO.Lm, 12. 38. Aw, 12.39.-x, 13.40.Gp 



The gauge-invariant, nonlocal, dynamical quark (GND) model [l| is one of a phenomenological non-local chiral quark 
r*^ I model that can be derived from QCD first principles by a series approximation [^. It appears as an approximate 
Oh' description for the cfFcctivc interactions among light quarks and pseudo-scalar mesons induced from the underlying 
^ ^ QCD after integrating out gluon and heavy-quark fields. In the low-energy region, the validity of the GND model is 

^ . tested by its resulting low-energy constants [l|-3 for the well-known Gasser-Leutwylcr chiral Lagrangian Q , which well 
f^ ' match the existing experimental data for pseudo-scalar meson physics. Considering that in obtaining the GND model 
the low-energy expansion is not taken, we expect this model has a larger range of application than the traditional low- 
energy region and are interested in its momentum-dependent behavior beyond the conventional low-energy expansion, 
►>. ' which reflects the interplay between low and high energies. In fact, the typical feature of the GND model is its 
OO ■ quark self-energy (or momentum-dependent quark mass) T,(—p'^) which represents the original idea of dynamical 
OO ' perturbation theory proposed by Pagels and Stokar [6|. AH dynamics, especially the momentum-dependent effects of 
OO , the model, are in the main effectively described by this quark self-energy. 

^^ ' To investigate the momentum-dependent effects in terms of a nonlocal model in general, the first obvious problem 

t:J" , one faces is its possible violation of the Ward-Takahashi identities (WTIs), as this happens for most nonlocal chiral 

^D ' quark models. In the literature, one way to solve the problem is to artificially revise the vertex [7|, |8|; another way 

^^ ' is to modify the model itself. The GND model belongs to this second approach, where the model is constructed in 

such a way that it is invariant under local chiral symmetry transformations for external current sources of light-quark 

fields. Except for the GND model, an earlier GNC model Q took the same tack but with a complex face factor 

introduced in the model. Because the local chiral symmetry is inserted into the model at inception, we expect the 

j^ WTIs to hold as a result. This expected validity of the WTIs was claimed but not explicitly shown in the GNC 

j^ , model, and not even mentioned in the original GND model. It is the principle aim of this paper to explicitly show 

that the GND model does satisfy the WTIs at the chiral limit. A by-product of this demonstration is that we not 

only give vector and axial-vector vertices, but also scalar, pseudo-scalar, and tensor vertices. In fact, due to their 

non-perturbative natures, these fundamental vertices have not been very well reported in the literature. The latest 

ansatz available for the vector vertex in QED is given in Refs.[l(j| which can be traced from the early Refs.Q. For 

vector and axial- vector vertices, the WTIs only constrain their longitudinal parts and cannot fix the transverse parts. 

For the remaining scalar, pseudo-scalar and tensor vertices, we even have no corresponding WTIs. Considering the 

vertices we obtained are from the unique GND model action and are constrained by inherent local chiral symmetry 

of the model, all different types of resulting vertices are at the same level as the approximate description of the 
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corresponding QCD ones. We will show that the vector and axial- vector vertices we obtained in the chiral limit are 
just the simplest versions satisfying the WTIs and have exactly the same form of those assumed by Ref.[lH; we call 
these two vertices and related scalar, pseudo-scalar, and tensor vertices the minimal Ward-Takahashi vertices. 

With proof of validity of the WTIs for the GND model and the resulting minimal vertices as our starting point 
to investigate the momentum-dependent effects, we take the light-cone pion distribution amplitudes (PDAs) as our 
next target of study in this paper. The reason these are chosen for discussion is that, at present, there is still no 
definite conclusion on whether PDA is in asymptotic-like form \12l, in Chernyak-Zhitnitsky (CZ)-like form |13, or in 
a flat-like form jl4| . With progress from experiments, ever more information and constraints have emerged on the 
PDAs providing good crosschecks between theoretical estimations and experimental data [15| . 

Theoretically, the model computation of PDAs is mainly through chiral quark models. Earlier calculations based 
on the local chiral quark model or local Nambu Jona-Lasino (NJL)-like models (see Ref.[lJ| and references therein) 
yielded the typical result that the lowest twist-2 PDA is in flat-like form in the chiral limit. In local chiral quark 
models, quarks have constant constituent masses and there is ultraviolet divergence in the resultant PDAs. To avoid 
the divergence, various regularization schemes are exploited which often yield confusing results. An improvement is 
to change to nonlocal chiral quark models (see Ref. J16l| and references therein), where a momentum-dependent quark 
mass or quark self-energy is arranged in such a way that it provides the model with a natural soft ultra-violet cutoff 
and results in finite PDAs. Some researchers believe that the momentum dependence of the quark self-energy will 
lead to twist-2 PDA deviating from flat-like form and generate correct end-point behavior [l7|. As we mentioned 
previously, a nonlocal chiral quark model usually violates the WTIs and the vertices are revised to avoid the defect. 
Now, our GND model automatically satisfies the WTIs, so there is no need to make such artificial modification of the 
vertices. 

Another technical problem met in nonlocal chiral quark model is that momentum integrations involving quark self- 
energies are usually difficult to complete, because usually these have on-shell external momenta which are time-like 
whereas conventionally loop momenta is space-like because the Wick rotation of the integration momenta occurs in 
Euclidean space. This mixed appearance of different kinds of momentum variables might cause variable p^ in the 
quark self-energy I](— p^) in the time-like re gion or even with imaginary cornponents. If there exists an analytical 
expression of I](— p^), such as for instantons 18| or just some simple ansatz [ll|, there will be no problem, because 
the extension to the time-like region or imaginary region is explicit. If I](— p^) however is defined as the solution of 
the Schwinger-Dyson equation (SDE), we meet a difficulty. The conventional numerical solution of the SDE is only 
defined in a space-like region; beyond that, the solution has still not been investigated well [l9[. In this paper, chiral 
perturbation is exploited to overcome this difficulty, with the finding that this perturbation works well. 

In the following, we first give a short review of the GND model, next compute vertices, proving that our vertices 
satisfy the standard WTIs, and finally discuss PDAs. We start with the generating functional of QCD, 

where ip and ip are light-quark fields (only u and d quark are taken as light quarks in this paper), ^ and ^ are 
heavy-quark fields with mass M, A" is the gluon field and G" its field strength. / and / are external sources for 
light-quark fields; J is the external current source for bilinear local light-quark fields. According to its 7 matrix 
structure, it can be decomposed into vector, axial- vector, scalar, pseudo-scalar, and tensor parts 

J(x) = fix) + ^{x)j5 - s{x) + ip{x)j5 + t'"'{x)a^^ . (2) 

The current quark mass m for light quarks are absorbed into the scalar source s{x); i.e. on the vacuum, all external 
sources vanish, except s{x) = m. After integrating out gluon and heavy-quark fields and integrating in the pseudo- 
scalar meson field, according to the discussion of Ref. [2| , the generating functional ([T]) can be approximated by 



SoNDiipn^i'n, Jn, In, In] = / ^^ 



/■^2, 



(4) 



4'i^[i9 + Jn - S(V )]Va2 + Inipn + Inipn 

i^n = [n'^PR + nPL\^ ^^-.^^[n^PR + nPL] In = [^^ Pl + nPR]i 7n=l[n^PL + nPii] (5) 

Jn = [n^PL + nPn] [J + i$\ [n^PL + nPii] = ^52 + ^j275 - sn + ipnlz + a^ut'^" (6) 

U = n' V^ = d^- iv^-^ , (7) 

where U and fi are unimodular pseudo-scalar meson fields, ■^ojV'o ^-nd In, In are rotated light-quark fields and 
corresponding sources, respectively, Jo is the rotated source for the rotated bilinear quark currents. I](— p^) is the 



self-energy for the rotated light-quark fields which satisfy the corresponding SDE. SGNu[''Pn,4'nyJniInjIn] is the 
action of the GND model; in the original paper [l|, / and / were not introduced, but an extra normalization term 
iTi\n[i0 + Jq] appeared. In Ref.[3|, this extra term was later proven to drop out. Integrating out the rotated 
light-quark and pseudo-scalar meson fields of ([3]), we obtain 



W[I,I,J] = — iTrhiL* ^ — / d'^xd'^yIn{x)D{x,y)In{y) -f pseudo-scalar meson loop corrections 



(8) 
(9) 



dn{x) 



for / 



0. Considering 



D~\x, y) = [z4 + Jo - S(Vj](5(a; - y) , 

where the pseudo-scalar field J7 in ([5]) must satisfy the stationary equation 

that in the low-energy region, we have shown in [^ that the Lagrangian of ([8]) just presumes the standard Gasser- 
Leutwyler chiral Lagrangian [5|, the stationary equation then can be replaced with the equation of motion (EOM) 
derived from the Gasser-Leutwyler chiral Lagrangian. For convenience in computation, we use below this chiral 
Lagrangian-induced EOM for 51, because it is more common and relatively simple. 



Neglecting pseudo-scalar meson loop corrections, the vertices in the chiral limit arc defined by the following 3-point 
Green's function 



[n\x)PL + n{x)Pi,]D{x,y)[n^{y)PL + n{y)PB] 



EE d'x'd^y'D{x,x')rnx',y',z)D{y'y) (10) 



SJ-iz) 
which yields 



,7=0,/=/=0 



Ttix,y,z) 






(11) 



j=o,/=/=o 



where i ^ S, P, V, A, T, t"" are the Pauli matrices in isospin space with r" = 1, and 

rs^^s^ rp^p^ J^^vl J%^al J^=il, 75 = 1 7P=»75 7y = 7^ 7^=7^75 It ^ ^r^" 

with the help of the following EOM results in the chiral limit 
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J=0, 1=1=0 



T''{l-5ao)+SaO 



(12) 
-^5{x-y), (13) 



where Bq — —^ {'fpip) /Fq is the p^-order low-energy constant of the Gasser-Leutwyler chiral Lagrangian related to the 
ratio of the quark condensate {tp'4') ^^'^ the square of the pion decay constant, Fq. With the help of ((T^ and (|13p . 
we can compute a series derivation of rotated sources Sii,p^,vn,aQ,i^ to un-rotatcd sources s,p,v,a,t. With these 
relations, we obtain the chiral-limit result 



rUx,y,z) = T^Six-z)d{y-z) 



^T.ni,ix, y, z) = r°cr^^(5(x - z)5{y - z) 



T1{x, y,^)=J ^^e-'9-+^Py+^(9-P)-f ^(p, q) t^PV^A 



(14) 
(15) 



Whereas the resulting scalar and tensor vertices are trivial, the pseudo-scalar, vector, and axial-vector vertices are 
non-trivial, their momentum space expressions being 



f?,(p,g)=zr° 



I -Boil -Sao) 



S](-g^) + S](-p^ 

(g-p)2 

7.-%^[S(-g^)-S(V)] 
q^ — p"^ 



75 



(16) 
(17) 
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p(V) + S(-p2)] 



75 • 



(18) 



Equations P7)) and P^ have exactly the same forms as those assumed by Ref. [Tl[ . From (|17p , the WTI for the vector 
vertex is 



{q-vYn 



V,iJ. 



^[S-\q)-S-\p)] 



S"\p)=f-n-P^ 



(19) 



where Sij)) is the Hght-quark propagator in momentum space with D{x,y) 
WTI for the axial- vector vertex is 



j=o 



/ ^e-'P-^S-^j). From ^, the 



{q-PnXu-^''''iS-\qh5+^5S-\p)] 



(20) 



The above two WTIs are standard and yield the minimal solutions for dTT]) and dUl). We call dUl), ^, ^ and ^ 
the minimal WT vertices and the corresponding GND model is the model which correctly generates these minimal 
WT vertices. In other words, if one prefers to trace the minimal WT vertices back to a source effective action, it is 
just the GND model action dH). 

Considering that the present GND model is only an approximation of the underlying QCD, we will in the future 
discuss its corrections in QCD. That will mean amending the above minimal WT vertices. 

We shall now discuss the light-cone PDAs, their definitions being as in Ref.|20| 



(p1V'(x)V(0)|0) = -^ / due'(i-")P-^ $(u,p,a;) , 



<P{u,p,x) 



y50(w)- ^'^5[</'p(") 



c^LuP^x"-^ — ] -I- high twist corrections 



p^— -1/ P^+?71^ 



(21) 
(22) 



where {p\ is a pion state with momentum p, ttIt^, rriu, rud are masses of the pion, u quark, and d quark respectively. For 
simplicity, we ignore the mass difference between u and d quarks setting them to the same current mass m = rUu = rnd- 
(j>{u) is the leading twist or twist-2 PDA, tpp{u) and 4>a{u) are sub-leading or twist-3 PDAs that correspond to the 
pseudo-scalar and the pseudo-tensor structures respectively. In our model, the operator ■i/;(a;)?/;(0) appeared in (j2ip is 
related to - i, ,.';„. , and its contribution to the pion matrix element relies on the part proportional to the pion field 
in the result. Parameterizing the pion field 11 by i7 = e*'^/^, the rotated source becomes 



Jn 



= [1 + *?75][-m + i$][l + 1^75] = -m - (zmH + ^^0)75 + 0{lY) 



J=-m, 1=1=0 



(23) 



Hence, 



{p\iJix)ij{om = (p\ 
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1 



*(pin(o)|o) 
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(2^' 



' - m - S((92) - [imlV + ^^0)75 
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- [m + S(-(72)]2 



756*^-^- 



(a;,0)[l + -n(x)75]|0) 

1 t 1 

^+7n-|-E(— g2) 2 f-(^-ni-T,[-{p-qY] 



,(24) 



where the first term, after expanding e*'''^ in terms of powers oi q- x and ignoring high twist 0{x'^) terms, contributes 
to the twist-3 pseudo-scalar PDA 4'p{u) at the end-point w = 0, whereas the second term can be changed to the form 
of (|2ip by the standard Feynman parametrization method. 
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{2^Y 
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r75,(25) 



where u is the standard Feynman parameter appearing in the well-known Feynman parameter integration and T,(q) 
is the abbreviated expression for m + J^l—q^). Conventionally, the Feynman parameter integration is used to deal 
with loop-momentum integration with constant mass in the perturbation theory; here, applied to the momentum 
integration involving momentum-dependent quark self-energy, it still works. One can easily cheek that if we take the 
quark self-energy appearing in p5|) back to constant mass as in the traditional NJL model, p5|) just becomes the 
standard Feynman parameterization formula. Indeed, this identification of the Feynman parameter with the PDA 
variable u not only endows the traditional mathematical Feynman parameter with a physical meaning, but also is 
valid in any kind of chiral quark or NJL-like models, as long as we have momentum integration of form (j25|) with two 
quark propagators (the form of vertex between two propagators is not important for this issue) , no matter whether 
the models are local or non-local. This result is independent of the GND model investigated in this paper; the GND 
model here is taken as an example to exhibit details of the computation . 

With the above Feynman parameterization, and after lengthy computations, we finally obtain $(it,p, x) introduced 
in (ini) as 



^$(u,p,a;) = t{p\UmO){S{u-0+) 



m 
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d^q 






3p2 



ih 



{2ttY {g2 + u(l _ y)p2 _ f;2(yp _ ^) _^ u[T?{up - q) - f^^ {q + p - Up)]}^ 

(p-q) 



75 



-5{u-\-) + 5{u-0+) + 



d'^q 



3p2 



h+m-i 



3 ^ 3p2 y p2 



(27r)4 {q2 ^ ^(1 _ y)p2 _ Y?{up - q) + u[S2(up - q) ~ Y?{q + p- up)]} 



-high twist terms 



where p^ — m^ and 



-J5 



(26) 



h 



phP + :^ (1 - 2u)]5^ + mE(Mp - g) - mS(g +p~ up)] + -p^[2(l - 2u)m + T,{up -q)-t{q+p- up)] ,(27) 



h = -x^py[-1^',Y][l^^{up- q) - m + -t{q + p - up)] , 

h = ("^ -f)q^ + [(1 - u)p-^{q+p- up)]{m+ -)[-uf+j:{up - q)] 

h = -t'iup- q) + {I - 2u)m--t,{q+p- up) . 



(28) 
(29) 
(30) 



In ()26p . we have dropped terms proportional to ^, since they belong to twist-4 terms. The delta function terms are 
just end-point terms that are from non-exponential p-x terms appearing in the original momentum integration ([25]) 
when we expand e"''^' in terms of powers oi q ■ x and apply Lorcntz invariance in decomposing its Lorentz structure. 
Since (j26p already has structure of (|22p . through comparison between the two equations, we can easily read out 
general expressions of PDAs 0(u), ippiu) and (/)cr(u) in terms of quark self energy. One can check that except the 
pure end point term in the first line of (j26p . all other terms satisfying symmetry u ^ 1 — u, since under combined 
transformation of u -i-> 1 — u and q O —q, the integrand of momentum integration is even for the second line and odd 
for the third line, respectively. This implies the corresponding symmetry of u -s-)- 1 — u for result PDAs. In the chiral 
limit, m = and p^ = mj = 0, we find ([26]) simplifies to 
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4(pin(o)|o) 
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2m '^ '^^^ 6 

where the four coefficients (t>5o,4'0:4'p,sa and ^g-.o are expressed as Euclidean-space momentum integrations 
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(34) 



Comparing the above result with ((22)) . we find that with the exception of end-point values for (j){u) symmetrically at 
u=0 and u=l and 4>p{u) non-symmetrically at u=0, PDAs are all independent of m {(j)p{u) even vanishes) and therefore 
take flat-like forms. Graphically, the pictures are that two symmetric infinitesimal narrow pillars at two end points 
appear above the flat-like form backgrounds for (^(w), one pure non-symmetric pillar at m = for 4>p{u) with zero 
background, and no pillar for (j)a-{u) with just flat-like form background. This result is completely due to the Feynman 
parameter description of PDAs, and therefore is valid for any type of chiral quark or NJL-like models, whether 
local or non-local (one can check by replacing quark self energy with constant mass that the end-point terms from 
(5— function for 0(w) and ippiu) exist even in local situations). This result is differs from those previously obtained in 
the literature, where researchers believed that the momentum dependence of the quark self-energy would force PDAs 
to deviate from flat-like forms in the chiral limit. Considering our result is analytical that does not rely on technical 
details of numerical computations and is valid for a large class of chiral quark models, we believe it is reliable. This 
result implies that, at least in the chiral limit, the momentum-dependent behavior of the quark self-energy, or more 
fundamentally the related non-locality of its interaction, causes no difference in the form of PDAs with constant mass 
or local interaction such as the NJL model. 

For the chiral limit result pip . (j){u) and (t>a{u) are ultraviolet divergent. The reason causing this divergence is 
that the quark self-energy of the GND model is for the rotated quark field. If the quark self-energy instead is for 
the un-rotated quark fields, as the conventional instanton model does, we need to replace the current quark mass m 
appearing in the numerator of the l.h.s. term of ([25]) with the quark self-energy and ignore the ^ term in the same 
numerator. This then will suppress the ultraviolet behavior of the momentum integration decreasing the ultraviolet 
divergences of PDAs . 

If we want to go beyond the chiral limit, the momentum integrations in the second and third line of (|26p are not 
easy to achieve. To finish the momentum integrations, we are used to rotate the momentum integration variable q into 
Euclidean space, although the external momentum p must be kept on the pion mass shell, p^ = m^. This will create 
some imaginary components; for example Ti[— (q — up)'^] = Y][—q'^ — u'^p'^ + 2uq-p] will become Ti[q'^ — u'^p'^ + 2u{iq'^p'^ — 
qE ■ p)] after a Wick rotation for integration variable g'' = (9°, q) — > {iq%^ qs)- Because the SDE now cannot provide 
us with a reliable quark self-energy beyond the space-like momentum region, we are not able then to directly compute 
the momentum integration in (j26p . One way to avoid this difficulty is to go back to the original constant quark mass 
case (or equivalently the NJL model situation) , where all momentum integration can be analytically finished except we 
need some momentum cutoff to regularize the integration. An alternative approach is to set an analytical expression 
for I](— (7^), as for instantons [l8| or introduce some simple ansatz 11|, then the momentum integration can still be 
finished, at least at the level of numerical computations. Considering that for the former information of the momentum 
dependence of the quark self-energy is lost, and the latter is constrained by specific choices of quark self-energy which 
also might not precisely describe its QCD behavior, we propose in this paper to expand the integrand in terms of 



powers of p 



and m. The underlying basis for this expansion is that when we go to higher-order terms in the 



expansion, according to (^5)) . typically we encounter a factor of [p-ql{q^ + S^)]^ or p^ /{q^ + T?) or mp-q/{q'^ + S^ 
In units of the QCD scale parameter Aqcd, the largest contribution comes mainly from the platform region of the 
quark self-energy in which E/Aqcd ^ 2 and (Z_e/Aqcd < 1 because quark self-energies above Aqcd decrease with 
momentum as 1/q^ [21 and contribute little. Considering that in our case Aqcd ^ 440MeV is much larger than the 
pion mass, i.e. p/Aqqd ^ itT'tt/^qcb ~ 1/3, and m < lOMeV, then the first factor [p-q/{q^ + Y?)Y ^ 10"^, the second 



factor p7(g^ + E^ 



10" 



and the third factor rap-q/{q^ + E^) ^ 10 ^] that is, each factor is at least an order of 



magnitude small. We expect chiral expansion will works well. With this analysis, the result of the detail computation 
gives up to first order. 
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We see that the first-order corrections include three parts 



(40) 
(41) 
(42) 
(43) 
(44) 
(45) 



• Corrections to chiral limit result which include: heights to two symmetric infinitesimal narrow pillars at two 
end points from (psi and to the flat-like form backgrounds from (piQ for (/)(u); heights to non-symmetric pillar at 
M = from (/)p.5io and to zero backgrounds from (ppi for (j)p{u); and height for flat-like form background from 
</>o-,io for 0<t(u). 

• Two symmetric infinitesimal narrow negative pillars at two end points from (fipju for 4>p{u) appear. 

• Asymptotic-like form concave type corrections [1^ proportional to u{l — u) from (pn and (/io-.ii to (j){u) and 
(j)a{u) respectively. 

To compare with literature's CZ-like form results, traditional double-hump structure of PDAs now for 0(w) is squeezed 
to two symmetric infinitesimal narrow pillars at two end points plus concave type asymptotic- like form above flat- 
like form background; for 4'p{u) is squeezed to two non-symmetric infinitesimal narrow pillars at two end points 
above a flat-like form background; for (j)a{u) is changed to concave type asymptotic- like form above fiat-like form 
background. We expect more complex non-asymptotic-like form corrections will show up in more higher order of our 
chiral perturbation expansion. 

If we further consider the normalization condition for (j){u) , 



du (j){u) 



then our results (PTj) . ([55]) and (pS)) imply that 



4(pln(O)|O)(205o + 2,^51 + 00 + <^io + TT-^ii) = 1 



6 



(46) 



(47) 



To obtain numerical results, we solve the SDE obtaining a quark self-energy as in Ref.[^ with model A and Aqcd = 

2 
440MeV, TOtt = 139.6MeV and m = ^^ = 12.3MeV, substitute the resulting quark self-energy into the above formulae 

for the coefficients, and finally perform the numerical computations. Considering that some momentum integrations 

are divergent, we take A = IGeV as the cutoff parameter in the Euclidean space momentum integration. Expressing 

all results in units of m,r, values for the coefficients are listed in Table|T] 

We see that the first-order corrections are orders of magnitude smaller than the leading order result. This verifies our 

conjecture that chiral perturbation will work well. 

To summarize, we have shown the GND model is a model which can correctly generate the minimal WT vertices. 
For PDAs. GND's result is similar to that of the NJL-like models, i.e. in the chiral limit except some end point 
pillars, PDAs take flat-like forms and non-flat effects of PDAs are due to nonzero pion and light-quark current mass 
corrections. We have shown that the variable u in PDAs is just the Feynman parameter of loop calculations of 



TABLE I: The obtained coefficients in unit of 10 



4>so 


(psi 
ml 


00 

ml 


010 

ml 


011 
ml 


4>v.so 
ml 


0p,«lO 

ml 


0p,511 

m3 


0pi 
m3 


4>cr,0 

ml 


0<T,1O 

ml 


0O-,11 


ml 


m2 


1.85 


-0.15 


10.13 


0.98 


-2.75 


82.22 


5.49 


-0.16 


0.86 


8.13 


-0.72 


-0.17 



standard perturbation theory and chiral perturbation works well for PDAs enabling quantitative estimates of PDA 
values. These results are valid not only for the GND model, but also for the larger class of chiral quark and NJL-like 
models. Further, the leading order nonzero pion and current quark masses corrections lead concave type asymptotic- 
like form modifications for 4>{u) with end-point pillars and (t>a{u) above the flat-like form backgrounds as a substitution 
of original double- hump structure of CZ-like form of PDAs. To force matching our results with other phenomenological 
ones especially the end-point behaviors, one can take present GND model chiral limit result as a starting point, ignore 
present corrections from pion and current quark masses, perform the QCD Gegenbauer evolution for PDAs as done 
in Ref.fll. 
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